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Many-electron dynamics induced by a symmetric monocycle electric-field pulse of large am¬ 
plitude is theoretically investigated in one- and two-dimensional half-filled extended Hubbard 
models on regular lattiees (i.e., without dimerization) using the exact diagonalization method 
for small systems and the Hartree-Fock approximation for large systems. The formation of a 
negative-temperature state and the change from repulsive interaetions to effective attractive 
interaetions are shown to be realized for a wide region of the field amplitude and the exeitation 
energy. For a nonnegligible intersite repulsive interaetion, the numerieal results are consistent 
with the fact that the phase separation between charge-rich and charge-poor regions is eaused 
by the corresponding effective attraetion. 


1. Introduction 

Physical phenomena eaused by intense optical pulses have attraeted much attention. For 
relatively weak pulses, various photoinduced phase transitions are known to take place in 
solids.Here, however, something qualitatively different is expected. Among the states that 
have been proposed for intense pulses, we focus on negative-temperature states.For sys¬ 
tems whose energy speetra are bounded above, negative-temperature states are realized when 
high-energy levels are oeeupied in sueh a way that the energy of the system is higher than that 
of an equilibrium state at infinite temperature.^ ’ To form such a state, we need to inerease 
the energy of the system (i.e., occupy high-energy levels selectively), suppressing the rise in 
the entropy (i.e., not merely occupying various energy levels). 

Two-level systems are typical ones for whieh negative-temperature states are discussed. 
Many moleeular materials have dimerized struetures in whieh bonding and antibonding or¬ 
bitals exist. They ean be regarded as interaeting two-level systems that are linked through in¬ 
terdimer eleetron transfers. The quasi-two-dimensional metal eomplex Et2Me2Sb[Pd(dmit)2]2 
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(dmit = l,3-dithiol-2-thione-4,5-dithiolate) is one of those strongly dimerized moleeular ma¬ 
terials for which photoinduced phase transitions have been experimentally* and theoreti- 
cally^^^^ studied. It has a charge-ordered ground state in which charge-rich and charge- 
poor dimers are regularly arrayed. The model used for this material is greatly simplified, 
keeping the dimerized structure, into a one-dimensional three-quarter-filled strongly dimer¬ 
ized extended Hubbard model with electron-phonon couplings, for which the formation of 
a negative-temperature state and the change from repulsive interactions to effective attrac¬ 
tive interactions are demonstrated in the band-insulator phase using the exact diagonalization 
method. This model is close to interacting two-level systems in that the intradimer transfer 
integral is much larger than the interdimer transfer integral. 

For regular systems without dimerization, Tsuji et al. have shown such changes of interac¬ 
tions in negative-temperature states for continuous waves,* half-cycle pulses, and asymmetric 
monocycle pulses"^ using the dynamical mean-field theory for the half-filled Hubbard model. 
For continuous waves, the modulation of the effective transfer integral is known in the context 
of dynamical localization^*^^* and given by the time average of the transfer integral with the 
Peierls phase factor. For very short pulses, the modulation is explained through the dynamical 
phase shift on the basis of the sudden approximation."^ 

The dynamical phase shift is proportional to the time integral of the electric field. In the 
sudden approximation, when the dynamical phase shift is n, the sign of the Peierls phase fac¬ 
tor is instantaneously inverted and the band structure is consequently inverted. For symmetric 
monocycle pulses, the time integral of the electric field is zero, so that the dynamical phase 
shift vanishes. For finite pulse widths, symmetric monocycle pulses make little contribution, 
and their effects are not sufficiently large to produce a negative-temperature state, in the dy¬ 
namical mean-field theory as described in Ref. 4. However, negative-temperature states are 
realized by symmetric monocycle pulses, as described in Ref. 12, and explained by the total- 
energy increments.The modulation of the effective transfer integral is defined through the 
total-energy increment and described by the zeroth-order Bessel function, as in the dynamical 
localization for continuous waves, in the two-level system. 

As mentioned above, the strongly dimerized model used in Ref. 12 is close to inter¬ 
acting two-level systems. This closeness may appear to be advantageous for the system to 
increase the energy suppressing the rise in the entropy. Regarding the fact that the repulsion- 
to-attraction conversion is hardly realized by symmetric monocycle pulses applied to the 
half-filled Hubbard model, in the dynamical mean-field theory as described in Ref. 4, is the 
distance of the regular Hubbard model from interacting two-level systems the reason? In 
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this paper, we show that the answer is no: the formation of a negative-temperature state and 
the change from repulsive interactions to effective attractive interactions are achieved by a 
symmetric monocycle pulse applied to the one- and two-dimensional half-filled extended 
Hubbard models on regular lattices, i.e., without dimerization. However, the total-energy in¬ 
crement is no longer described by the zeroth-order Bessel function. 

Photoinduced superconductivity is discussed in Ref. 3 using the Hubbard model that has 
only the on-site interaction U. Real substances always have intersite interactions, which are 
important for the exciton effect. From the viewpoint of the possible inversion of interactions, 
the intersite interactions are important because they can cause phase separation once they are 
converted into effective attractive interactions. In this paper, we consider extended Hubbard 
models, which have the nearest-neighbor interaction V in addition. 

For phase separation, many theoretical studies have been performed using the one- 
dimensional^^ and two-dimensional^’ t-J models as well as the one-dimensional,^^’’ two- 
dimensional,’^ and even higher-dimensional’"^ extended Hubbard models at^* and near’"^ half 
filling, at quarter filling,and at general fillings.Many of them have been motivated by 
a desire to clarify the mechanism of high-temperature superconductivity. A superconducting 
state is indeed anticipated near the instability toward phase separation. It is now known that 
phase separation is realized in a rather wide region in the plane spanned by the interaction 
strength and the filling. In general, when the intersite interaction is attractive, phase sepa¬ 
ration is easily achieved. Thus, when the system has an intersite interaction and the 
interactions are expected to be inverted by an intense pulse, one must consider the possibility 
for phase separation. In this paper, we employ the exact diagonalization method for small 
systems and the Hartree-Fock approximation for large systems to show that phase separation 
is indeed possibly induced by an intense symmetric monocycle optical pulse. 

2. Extended Hubbard Model on Regular Lattice 

Many-electron dynamics after strong photoexcitation is studied in the half-filled extended 
Hubbard model, 

i ^ ' ' (ij) 

where creates an electron with spin cr at site i, ni^cr = and n,- = 2o- ^Ua-- The 

transfer integral is denoted by to- The parameter U represents the on-site repulsion strength, 
and V represents the nearest-neighbor repulsion strength. For the one-dimensional lattice, we 
use the 12-site periodic chain. For the two-dimensional lattice, we use the 4 x 3-site system 
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Fig. 1. 4 X 3 lattice with periodic boundary condition. 


with the periodic boundary condition shown in Fig. 1. 

Before the photoexcitation, the system is in the ground state obtained by the exact diago- 
nalization method unless stated otherwise. Photoexcitation is introduced through the Peierls 
phase 




exp 


le 

-^rij ■ A(t) 
He 
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where rij = rj - r, and r, being the location of the zth site. The right-hand side is substituted 
into Eq. (1). For symmetric monocycle electric-field pulses, we use the time-dependent vector 
potential 

A(t) = — [cos(mt) - 1] 6(t)9i— - n , (3) 

io \a) I 

where F = (F, F) is polarized in the (1,1) direction on the square lattice {F = F on the chain) 
and the central frequency co is chosen to be nearly resonant with the optical gap unless stated 
otherwise. The time-dependent Schrodinger equation is numerically solved by expanding the 
exponential evolution operator with a time slice dt=0.0\ to the 15th order and by checking 
the conservation of the norm.^^ 


The kinetic energy is defined as the expectation value of the first term in Eq. (1). The 
total energy is the expectation value of Eq. (1), which becomes zero in equilibrium at infinite 
temperature. We set to=1 as a unit of energy. The interaction parameters U and V are varied. 
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The intersite distance is set to be a in any direction. The time-averaged double occupancy is 
calculated by 

1 nts+tw 

= — I WO I I , (4) 

Jts 

with ts = 5T, tw = 45r, and T = Injco being the period of the oscillating electric field. The 
other time-averaged quantities are calculated likewise. 

The systems to which the exact diagonalization method can be applied are quite small, 
and their nonequilibrium dynamics basically from the ground state can be discussed. On the 
other hand, the dynamical mean-field theory can be used to discuss the thermodynamic limit 
at finite temperatures from the beginning, but it uses a mapping to a single-impurity problem, 
so that events in the momentum space such as Fermi surface nesting or the Umklapp pro¬ 
cess are included indirectly through the self-energy that is in the second order with respect 
to electron transfers. Thus, it often underestimates the tendency toward an insulator. Here 
only, we use a half-cycle pulse and compare numerical results obtained by the present ex¬ 
act diagonalization method with those in the dynamical mean-field theory for the half-filled 
Hubbard model (V = 0).^ The details are shown in the Appendix. In both approaches, we 
find (i) that negative-temperature states appear when the Peierls phase shift after the photoex¬ 
citation is near an odd multiple of n and U /(4to) is not too large, and (ii) that, as the pulse 
width increases and U /(4to) increases, the region where negative-temperature states appear is 
shifted toward a larger Peierls phase shift. These findings suggest that, for the formation of a 
negative-temperature state, it does not matter whether the initial state is metallic as in Ref. 4 
or insulating as in the present study for U 

3. States after Symmetric Monocycle Pulse Excitation in Hubbard Model 

Because the results for the one- and two-dimensional cases are similar, we show them 
in the two-dimensional case. Time-averaged quantities are shown in Fig. 2 as functions of 
eaFjTiio for U =4, 6, and 8. Around eaFlfito=3J, the time-averaged kinetic energy is positive 
[Fig. 2(a)]. This indicates that the momentum distributions at -2to[cos(kx) + cos(ky)] > 0 are 
larger than those at -2to[cos(kx)+cos(ky)] < 0. Note that a symmetric monocycle pulse is used 
and consequently the Peierls phase returns to the initial value. This is not caused by the band- 
structure inversion through the dynamical phase shift; Umklapp scattering is considered to be 
responsible. This change in the momentum distributions is accompanied by a large increase 
in the total energy. The total energy here is the sum of the kinetic energy and U times the 
double occupancy relative to its noninteracting value of 0.25 at half filling. 
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Fig. 2. (Color online) (a) Time-averaged kinetic energy and (b) time-averaged double occupancy 
as functions of eaFITiu for different values of U. The other parameters are V -0,u - 2.4 for U = 4, cu = 3.43 
for U - 6, and a> - 4.82 for t/ = 8. 


The time-averaged double oeeupaney indeed beeomes larger than 0.25 for 

an eaFIfico region ineluding the region for positive time-averaged kinetie energies around 
eaF/fi(i)=3J [Fig. 2(b)]. It appears as if these nonequilibrium states possess an attraetive 
on-site interaetion. This is eaused by the faet that the energy supplied by photoexeitation is 
distributed to both the kinetie energy and the interaetion energy. When the kinetie energy 
is positive and > 0.25, the total energy is positive, i.e., it is larger than the value 

taken by an equilibrium state at infinite temperature. This total energy is treated as an order 
parameter in Ref. 4. 

In the band-insulator phase of the one-dimensional three-quarter-filled strongly dimerized 
extended Peierls-Hubbard model, the region of eaFIhco where negative-temperature states 
appear is almost independent of U, and in the negative-temperature state (i.e., its 

effeetive on-site attraetion) inereases with U Here, however, the situation is not so simple. 
The value of eaF/Hco where the loeal maximum of the time-averaged kinetie energy appears 
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Fig. 3. (Color online) (a) Regions where time-averaged kinetic energy is positive, denoted by “KE>0”, and 
(b) regions where {(n,,fn,,j,)) is larger than 0.25, denoted by “DO>0.25”, for f/=4 in plane spanned by eaFItioj 
and w/Fopt with Fopt being the optical gap. 


or where that of {{n^rii ])) appears depends on U. Furthermore, in the band-insulator phase of 
the model above, the energy flow into the eleetron-eleetron interaetion terms is quite small. 

In the present ease, the insulating ground state is eaused by the eleetron-eleetron interaetion 
term so that the energy flow into this term is relatively large. This is quantitatively diseussed 
later in Seet. 4. 

Sueh negative-temperature states aeeompanied by inverted on-site interaetions appear in 
a wide m region. Their appearanee is not limited to the ease where m is nearly resonant with 
the optieal gap. For U =4, the regions where the time-averaged kinetie energy is positive and 
the regions where {{n^rii^i)) is larger than 0.25 are shown in Figs. 3(a) and 3(b), respee- 
tively. Here, the plane is spanned by eaFIfioj and m/^opt with ^opt being the optieal gap. For 
eaeh m/^opt, the upper and lower bounds of the region(s) are denoted by inverted triangles 
and triangles, respeetively. For m=0.6£’opt and O-V^opt, eaFIfioj must be larger than 6 for the 
time-averaged kinetic energy to be positive, but ((n,,tn,,j,)) > 0.25 is widely achieved. For 
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V 


Fig. 4. (Color online) Double occupancy (n/jn/^i) and nearest-neighbor charge-density correlation (n,n,+i) as 
a function of V in ground states with t/ = 8 and -8. 


a»=1.4£’opt and LS^opt, eaFlho) must be less than 4 for a positive time-averaged kinetie en¬ 
ergy and for {{n^nii)) > 0.25. In this sense, the resonant eondition is optimum for these 
two phenomena and eonsequently for a negative-temperature state to be formed. The regions 
of eaFIfioj for > 0.25 eontain those for positive time-averaged kinetie energies. 

Therefore, the former regions are wider than the latter. In the one-dimensional ease, negative- 
temperature states aeeompanied by inverted on-site interaetions also appear in a wide U,oj 
region. 

4. States after Symmetric Monocycle Pulse Excitation in Extended Hubbard Model 

If both the on-site and intersite interaetions are inverted, we expeet that the effee- 
tive attraetive interaetions will eause phase separation into a region of high-eharge-density 
sites and a region of low-eharge-density sites, although their boundaries fluetuate quantum- 
meehanieally. How easily the phase separation is realized depends on the dimensionality 
through the nature of the fluetuating boundary. Thus, below we eompare the one-dimensional 
ease, where the boundary eonsists of two points, with the two-dimensional ease, where the 
boundary eonsists of a line. 

4.1 One-dimensional case 

Before the analysis of photoindueed states, we show the V dependenee of the double oe- 
eupaney the nearest-neighbor eharge-density eorrelation in the ground 

states with f/ = 8 and -8 in Fig. 4. When the interaetions are repulsive, the spin-density- 
wave (SDW) eorrelation is dominant for V < t//2 -b 6, and the eharge-density-wave (CDW) 
eorrelation is dominant for V > f//2 -b 6. Here, <5 is a small positive quantity originating 
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Fig. 5. (Color online) (a) Time-averaged kinetic energy, (b) time-averaged double occupancy , 4 ,)), 

(c) time-averaged nearest-neighbor charge-density correlation ((n,n,+i)), and (d) total energy as functions of 
eaFIho) for different values of V. The other parameters are t/ = 8 , cu = 4.93 for V = 0, o) - 4.86 for V = 1, 
o) - 4.56 for V - 2, and a> - 3.72 for V = 3. 


from spin fluctuations that approaches zero in the to 0 limit. When the SDW eorrelation is 
dominant, is small and 1 ) is almost 1. When the CDW eorrelation is dominant, 

(niHi+i) is small and j^) is about 0.5. When the interaetions are attraetive, phase separa¬ 
tion is antieipated. In this ease, (n,--pn, j^) is about 0.5, whieh is the average of (n, ^ 1 on 

the high-eharge-density sites and - 0 on the low-eharge-density sites. Meanwhile, 

is about 1.67, whieh is the average of {ntni+i) ^ 4 on the bonds between high-eharge- 
density sites and {nirii+i) ^ 0 on the other bonds. In Fig. 4, we show the results for V > -0.3 
beeause the states for V < -0.4 are numerieally unstable. 

Time-averaged quantities for U=8 with dilferent values of V are shown in Fig. 5 as fune- 
tions of eaFIhco. In all the eases shown here, there is a region around eaF/fi.(o=3 for positive 
time-averaged kinetie energies, although this region is narrowed with inereasing V [Fig. 5(a)]. 
In nearly the same eaFIfia> region, the time-averaged double oeeupaney {{n^n^)) beeomes 
larger than 0.25 [Fig. 5(b)]: the elfeetive on-site interaetion is attraetive. On the other hand. 
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where the quantities shown in Figs. 5(a) and 5(b) are at loeal minima around eaFlfico=A.5, the 
time-averaged nearest-neighbor eharge-density eorrelation ((n,n,+i)) is at a loeal maximum 
[Fig. 5(e)], and the total energy is at a loeal minimum [Fig. 5(d)]. The eorresponding time- 
averaged state approaehes the ground state, although the former is not so elose to the latter 
eompared with the eorresponding time-averaged state in the one-dimensional three-quarter- 
filled strongly dimerized extended Peierls-Hubbard model, whose total-energy inerement is 
deseribed reasonably well by the zeroth-order Bessel funetion multiplied by the sine fune- 
tion^^ as in the two-level system.'^ 

Where the quantities shown in Figs. 5(a) and 5(b) are at the maxima around eaFlfioj=3, 
((n,n,+i)) is at a loeal maximum (minimum) for V=2 and 3 (y=0 and 1) [Fig. 5(e)], and the 
total energy is positive and at the maximum [Fig. 5(d)]. The state realized here is a negative- 
temperature state. The faet that, for y=2 and 3, both ((n/jn,^)) and ((n;n,+i)) are at loeal 
maxima in the negative-temperature state indieates that phase separation is realized here. 
To see how roughly the negative-temperature state is interpreted as the ground state with the 
inverted interaetions, we eompare the eorrelations in the ground states for U =8 and y > 0 and 
for t/=-8 and y < 0 (Fig. 4). For small and positive V (y=0.1 and 0.2), is deereased 

by the inversion of the interaetions. This is in eontrast to the eases of y > 0.3, where both 
(ni -fn/ i) and («;«;+!) are inereased by the inversion. The loeal minimum of ((«,«,+ 1 )) around 
eaF/fico=3 for y=l [Fig. 5(e)] should be eaused by the faet that V is small and the inversion 
of the interaetions is ineomplete. As a matter of faet, ((n,n,+i)) at the loeal maximum is 
smaller (only slightly larger) than (niHi+i) in the ground state for y=2 (y=3). Thus, the overall 
behaviors of the eorrelations in the negative-temperature states are eonsistent with those in 
the ground states with ineompletely inverted interaetions. Compared with the inverted on¬ 
site interaetion, whieh makes the double oeeupaney larger than 0.25, the inverted intersite 
interaetion is less effeetive. Note that, in the ground state, the value of | y | required for phase 
separation (U, V < 0) is generally mueh smaller than that for the CDW ([/, y > 0). In the 
thermodynamie limit and in the to ^ 0 limit, the former approaehes zero but the latter is 
proportional to U with a dimension-dependent eoefheient. From this viewpoint, even if the 
inversion of the intersite interaetion is ineomplete, photoindueed phase separation is feasible. 

4.2 Two-dimensional case 

First, we show the V dependenee of the double oeeupaney and the average 

nearest-neighbor eharge-density eorrelation j{nini+^ nini+^), where x and y are the unit vee- 
tors along the x- and y-axes, respeetively, in the ground states with [7 = 8 and -8 in Fig. 6. 
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V 


Fig. 6. (Color online) Double occupancy {n^rii^i) and average nearest-neighbor charge-density correlation 
^{riitii+s, + nifii+y) as a function of V in ground states with U - 8 and -8. 


When the interaetions are repulsive, the SDW eorrelation is dominant for V < U/4 + 6, 
and the CDW eorrelation is dominant for V > U/4 + 6. Here, 6 is not so small beeause the 
present 4x3 lattiee is ineommensurate with the CDW with the wave veetor of {n,n). In the 
thermodynamie limit and in the to ^ 0 limit, 6 approaehes zero. Compared with the one¬ 
dimensional ease, the numerieal instability for attraetive interaetions is greatly suppressed. In 
the ease of phase separation, is about 0.5 and ^{nirii+s: + riiHi^^) is about 1.17, whieh 

is the average of j{nini+^ + nini+^) ^ 4 on the bonds between high-eharge-density sites and 
^ 0 on the other bonds. Note that the weight of the boundary bonds between 
the high-eharge-density sites and the low-eharge-density sites is mueh larger than that in the 
one-dimensional ease, so that the nearest-neighbor eharge-density eorrelation here is substan¬ 
tially smaller than that in the one-dimensional ease. In other words, the boundary degrees of 
freedom are mueh larger in the 4x3 lattiee than in the 12-site ehain, so that the tendeney 
toward phase separation is mueh weaker in the 4x3 lattiee. 

Time-averaged quantities for [7=8 with different values of V are shown in Fig. 7 as fune- 
tions of eaF/hto. Compared with the V/2 dependenee of the eorresponding quantities in the 
one-dimensional ease, the present V dependenee is generally small, whieh is eonsistent with 
the faet that the tendeney toward phase separation is weaker in the two-dimensional ease. 
Then, we show the region of 1 < eaF/Uco < 4 to enlarge the view around the loeal maxima of 
the total energy. The eaF/fio) region for positive time-averaged kinetie energies is shifted to 
the smaller eaF/fio) side with inereasing V [Fig. 7 (a)]. In a relatively wide region ineluding 
the eaF/fioj region above, ((n,,tR,,j,)) beeomes larger than 0.25 [Fig. 7(b)]: the effeetive on-site 
interaetion is attraetive whenever the time-averaged kinetie energy is positive. 
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Fig. 7. (Color online) (a) Time-averaged kinetic energy, (b) time-averaged double occupancy (c) 

time-averaged nearest-neighbor charge-density correlation + n,n/+y)), and (d) total energy as functions 

of eaFlhu for different values of V. The other parameters are t/ = 8, w = 4.82 for V - 0,a) - 4.75 for V = 0.5, 
Cl) = 4.64 for y = 1, and co - 4.47 for V - 1.5. 


Where the quantities shown in Figs. 7(a) and 7(b) are at the maxima around eaF/tico=3.5, 
^{{nirii+jc + nini+0) is at a local maximum (minimum) for y=l and 1.5 (V=0) [Fig. 7(c)], and 
the total energy is positive and at the maximum [Fig. 7(d)]. Negative-temperature states are 
realized here in a wide eaFlhco region because the U term assists the total energy to become 
positive, while the V term does not assist it at all (^((n;n,+jt -l-n;n,+j))) is always smaller than 1). 
The fact that, for y=l, 1.5 (shown here), 2, 2.5, and 3 (not shown here), both and 

^{{TiiTii+y- -b nini+0) are at local maxima in the negative-temperature state indicates that phase 
separation is realized around eaFl%a)=3.5 in the two-dimensional case. Then, we compare 
the correlations in the ground states for t/=8 and y > 0 and for U=-S and y < 0 (Fig. 6). 
For small and positive V (y=0.1 and 0.2), i((n;n,+jE 4- nini+^)) is decreased by the inversion 
of the interactions. The absence of an apparent local maximum of ^((n,n,+jE -l- n,n,+ 5 i)) around 
eaFlfiaj=3.5 for y=0.5 [Fig. 7(c)] should again be caused by the fact that V is small and 
the inversion of the interactions is incomplete. As a matter of fact, ^+jt + at the 
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local maximum is smaller (larger) than ^{nini+^ + nini+^) in the ground state for V=l and 1.5 
(V=2, 2.5, and 3) (not shown). Again, the overall behaviors of the correlations in the negative- 
temperature states are consistent with those in the ground states with incompletely inverted 
interactions. Even though f ((n,n,+^ -l- at the local maximum is still smaller than 1 in 

all cases in the present small system to which the exact diagonalization method is applicable, 
a strong photoexcitation may lead to phase separation in much larger systems. 

Recall that the formation of a negative-temperature state is generally obstructed by a 
large rise in the entropy. The photoinduced phase separation is not obstructed by such an 
effect. However, when we employ attractive interactions (U, V < 0) and strongly excite 
phase-separating states, no negative-temperature state is realized (not shown). The total en¬ 
ergy and the time-averaged kinetic energy remain negative for any eaF/fioj under both reso¬ 
nant and nonresonant conditions. The time-averaged double occupancy is always larger than 
0.25. These facts indicate that it is difficult to increase the total energy of an initially phase- 
separating state without a large rise in the entropy. The entropy is easily raised by exciting 
low-energy collective states corresponding to the fluctuating boundary. 

4.3 Hartree-Fock approximation for two-dimensional case 

It is difficult to treat larger systems at half filling for long periods of time by the exact diag¬ 
onalization method. However, in order to judge whether phase separation is really achieved, 
we need to study the photoinduced dynamics of larger systems anyway. Hereafter, we use 
the Hartree-Fock approximation for larger systems with a periodic boundary condition on the 
square lattice and compare the results for different system sizes. The ground state is an SDW 
state. The central frequency to is chosen to be nearly resonant with the optical gap in the time- 
dependent Hartree-Fock approximation. In expanding the exponential evolution operator, a 
time slice of dt = 10“^ is used. The time averaging is performed using Eq. (4) with ts = lOOT 
and F = 50r. At first, in order to compare results obtained by these methods, we show the 
same quantities as in Fig. 7 as functions of eaF/hto. Here, the kinetic and total energies are 
divided by the number of sites for the comparison of results with different system sizes. The 
data shown here are obtained for the 10x10 lattice, and they coincide with the results for the 
14x14 and 16x16 lattices within the symbol sizes. Although the value of eaF/hto at which a 
local maximum/minimum of each quantity appears deviates from the corresponding one ob¬ 
tained by the exact diagonalization method, these results obtained with the different methods 
show qualitatively similar eaF/hto and V dependences. 

In a region around eaF/hto=2.A, the time-averaged kinetic energy and the total energy 
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Fig. 8. (Color online) (a) Time-averaged kinetic energy per site, (b) time-averaged double occupancy 
(c) time-averaged nearest-neighbor charge-density correlation + n,n,+y)), and (d) total en¬ 

ergy per site as functions of eaFItiu obtained in Hartree-Fock approximation for different values of V. The other 
parameters are t/ = 8, w = 7.0 for V = 0, at - 6.6 for V — 0.5, a> - 6.2 for V = 1, and a> - 5.8 for V = 1.5. 


are positive [Figs. 8(a) and 8(d)]. However, their V dependenee is opposite to that obtained 
by the exact diagonalization method, although the V dependence itself is rather weak. The 
total energy shows a local maximum near eaF/fi.(o=2A, where becomes larger 

than 0.25 [Fig. 8(b)]. The effective on-site attraction increases with V. Around eaF/fito=2A, 
^{{TiiTii+sc + TiiTii^^)) is at a local maximum (minimum) for F=0.5, 1, and 1.5 (y=0) [Fig. 8(c)]. 
Now, ^{{riirii+i + nini+^)) exceeds 1, but its maximum value is slightly over 1. The behavior 
around eaFltia)=2A is consistent with the appearance of phase separation, but ^{{nini^jc + 
nini+0) is considerably smaller than that in the phase-separating ground state with the inverted 
on-site and intersite interactions. This fact will be clarified below when the space and time 
dependences are discussed. The maximum value of ^((n,n,+jf -b n,n,+y)) is insensitive to the 
system size, at least up to the 16x16 lattice. As far as the quantities shown here are concerned, 
they should already be close to those in the thermodynamic limit. However, the space and time 



















Hartree-Fock, U=8, V=^, eaF/ha)=2A 



Fig. 9. (Color online) Snapshots of charge density (n,) at (a) f=60 on lOx 10 lattice, (b) f=80 on lOx 10 lattice, 
(c) f=160 on 14x14 lattice, and (d) f=200 on 14x14 lattice obtained by Hartree-Fock approximation for t/ = 8, 
V - I, u - 6.2, and eaFlhco - 2.4. 


dependences can be sensitive to the system size. 

Snapshots of the charge-density distribution (n,) are shown in Fig. 9 at different times 
and on the 10x10 and 14x14 lattices. The contour lines for (n,)=0.7, 1.0, and 1.3 are shown 
on their bases. The initial condition before the photoexcitation is spatially uniform, so that 
{«;)=! at t=0. However, the translational symmetry becomes spontaneously broken after the 
photoexcitation. Then, the charge-density distribution (n,) fluctuates spatially and temporally. 
At first, a large number of small high-density regions and low-density regions appear. As time 
proceeds, some of them merge into larger high-density regions and low-density regions. On 
the 10x10 lattice, there are two high-density regions at t=60 [Fig. 9(a)], and they have already 
merged into one at t=80 [Fig. 9(b)]. On the 14x14 lattice, there are two high-density regions 
at 1=160 [Fig. 9(c)], and they are about to merge into one at 1=200 [Fig. 9(d)]. For the 16x16 
lattice, the time when the high-density regions finally merge into one is later than that for the 
14x14 lattice. 
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Such a size dependence of domain dynamics is actually reflected in a small but sig¬ 
nificant size dependence of ^-l- For t/=8, y=l, and eaF/fia>=2A, where 

“I” f'l shows a maximum, -i- — 1.079 on the 10x10 lattice 

where the high-density regions have already merged into one before ts (the time averaging), 
^((riini+x + n,n,+y))= 1.070 on the 14x14 lattice where the high-density regions have not yet 
merged into one at ts + F, and ^((n,n,+jE -l- n;n;+ 5 i))= 1.074 on the 16x16 lattice where the situ¬ 
ation is similar to that on the 14x14 lattice. On the 14x14 lattice, the existence of more than 
one domain makes the boundary contribution larger and the above correlation smaller than 
those on the 10x10 lattice. On the 16x16 lattice, the larger system size makes the bound¬ 
ary contribution smaller and the above correlation larger than those on the 14x14 lattice. The 
size dependence of the time-averaged correlation above thus originates from the timing of the 
merging and the relative weight of the boundary degrees of freedom. In any case, the time- 
averaged correlation above is almost saturated for these system sizes and the size dependence 
is small. 

The charge density in the high-density regions is not as high as (n,) - 2, and that in the 
low-density regions is not as low as (n,) - 0. This is regarded as being due to a rise in the 
entropy. The rise in the entropy is not so large as to obstruct the formation of a negative- 
temperature state, but it is significant, so that the charge disproportionation is small. This 
is the reason why the quantities shown in Fig. 8 are insensitive to the system size and the 
reason why the maximum value of -i- n,n,+ 5 i)) is slightly over 1. Even though the 

charge disproportionation is small. Fig. 9 shows a nonequilibrium phase-separating state, 
which is quantitatively different from the phase-separating ground state with the completely 
inverted on-site and intersite interactions. The inversion of the intersite interaction is incom¬ 
plete, namely, the effective intersite interaction becomes negative, but its magnitude is smaller 
than the initial value. Therefore, the previous results obtained by the exact diagonalization 
method for a small system are consistent with the present Hartree-Fock results and the antic¬ 
ipation of photoinduced phase separation in sufficiently large systems. 

5. Conclusions and Discussion 

Many-electron dynamics in the one- and two-dimensional half-filled extended Hubbard 
models after the application of a symmetric monocycle electric-field pulse are calculated us¬ 
ing the exact diagonalization method for small systems and the Hartree-Fock approximation 
for large systems. The appearance of negative-temperature states is demonstrated by their 
total energies higher than that in an equilibrium state at infinite temperature. They are char- 
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acterized by inverted on-site and nearest-neighbor interaetions, whieh are roughly estimated 
from the field-amplitude dependenee of time-averaged eorrelation funetions. The important 
point is that, even by applying symmetrie monoeyele pulses and even in models on regular lat- 
tiees without dimerization, negative-temperature states are realized in a wide region spanned 
by the field amplitude and the photoexeitation energy. Under nearly the same eonditions, the 
time-averaged double oeeupaney is larger than that of the noninteraeting equilibrium state, 
and thus the on-site interaetion is effeetively inverted. In addition, if the nearest-neighbor 
repulsion is not very small, the time-averaged nearest-neighbor eharge-density eorrelation 
shows a loeal maximum as a funetion of the field amplitude. Although the inversion of the 
interaetions is ineomplete, the intersite interaetion is also inverted. These field-amplitude and 
intersite-interaetion dependenees of the time-averaged eorrelation funetions are eonsistent 
with the antieipated photoindueed phase separation. 

The diffieulty of forming a negative-temperature state by applying a symmetric mono¬ 
cycle pulse in the dynamical mean-field theory for the half-filled Hubbard model"^ may be 
caused by the fact that Umklapp scattering is included indirectly through the self-energy in 
the theory, where the momentum space is basically irrelevant. The fact that the kinetic en¬ 
ergy becomes positive and consequently the momentum distribution does not return to the 
initial one after intense irradiation of a symmetric monocycle pulse suggests the importance 
of the Umklapp process. Thus, it would be desirable to treat the momentum space or the real 
space in a direct manner. The present approach is expected to be suitable for the discussion 
of negative-temperature states. Photoindueed superconductivity is proposed using the Hub¬ 
bard model.^ Real substances always have intersite repulsions, so that photoindueed phase 
separation can be realized once the intersite interactions are effectively inverted: not only two 
electrons but many electrons attract each other. As to phase separation, a spatially extended 
region must be treated to account for the gradual space dependence of the charge density. In 
this context, it would be better to supplement the exact diagonalization study with Hartree- 
Fock studies, as in this paper. Generally, a superconducting state with a high transition tem¬ 
perature is expected to be realized near some electronic instability. Phase separation has been 
discussed extensively in this context. Namely, in the vicinity of a phase-separating state, su¬ 
perconductivity with an electronic mechanism has been anticipated. As the possibility of 
photoindueed phase separation is enhanced, that of photoindueed superconductivity should 
also be enhanced. 

In our numerical calculations, as the system size increases, it requires a longer time for 
high-density or low-density regions to merge into one. This is due to the initial condition 
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in which the charge density is spatially uniform. In real substances at finite temperatures, 
electron-phonon interactions exist and thermal lattice fluctuations should exert a random force 
on electrons, so that high-density or low-density regions should easily appear, be activated to 
move around, and encounter similar-density regions that they merge with. Unless the system 
is biased externally, however, such regions can appear anywhere in the sample. We expect 
that an external bias such as a DC electric field would assist macroscopic phase separation 
with an anomalous dielectric property. Although the electron configuration in the ground state 
is determined so that the total energy is minimized, that in the photoinduced state with the 
largest total energy is determined in the opposite manner: the instability is a key because 
the realized state is the most unstable state among the photoinduced states. We expect that 
the negative-temperature state in real substances would decay rapidly through every chan¬ 
nel, interacting with phonons, photons, and electrons in the environment. The lifetime of the 
negative-temperature state would be short especially when the system has a complex structure 
and interacts with many degrees of freedom. 
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Appendix: Negative-Temperature States Induced by Half-Cycle Pulse 

Here only, we consider the half-cycle cosine pulse used in Ref. 4, which is in the present 
formulation written as 


fic 

A(t) = -(A, A) 

ea 


t 1 . I2nt\ 


9(t)e(T -1) 


(AT) 


T 2n 

with T being the pulse width. It is applied to the half-filled Hubbard model (V = 0) on the 
4x3 lattice. We use the notation A, as in Ref. 4, which is the magnitude of the Peierls phase 
shift for the electron transfer in any direction. 

The energy dependence of the density of states is different from that used in Ref. 4. 
Then, here, we take the half-bandwidth in the noninteracting case, 4to, as a unit of energy. 
The dimensionless repulsion strength is U /(4to) and the dimensionless pulse width is At^T. 
The region where the total energy is positive after the photoexcitation, i.e. where negative- 
temperature states appear, is surrounded by the dashed line for the delta-function pulse (r —> 
0) and by the dotted line for the pulse with width AtoT = 10 in Fig. AT. For the delta-function 
pulse, this region is symmetric with respect to the A = n line. For the pulse with finite width. 
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Fig. A4. (Color online) Region where total energy is positive, around point {U, A)=(0, n) in plane spanned 
by t//(4fo) and A, for Hubbard model driven by half-cycle cosine pulse with width r —» 0 or 4foT =10. 


this region is shifted to the larger A side as [//(4fo) inereases. These results are eonsistent 
with those in the dynamieal mean-field theory."* A negative-temperature state is formed by 
the Peierls phase shift as well in the present system. 
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